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Abstract 

We consider the M/M/l-PS queue with processor sharing. We study 
the conditional sojourn time distribution of an arriving customer, con- 
ditioned on the number of other customers present. A new formula is 
obtained for the conditional sojourn time distribution, using a discrete 
Green's function. This is shown to be equivalent to some classic results 
of Pollaczeck and Vaulot from 1946. Then various asymptotic limits 
are studied, including large time and/or large number of customers 
present, and heavy traffic, where the arrival rate is only slightly less 
than the service rate. 
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1 Introduction 



One of the most interesting service disciplines in queueing theory is that of 
processor sharing (PS). Here every customer in the system gets an equal 
fraction of the server or processor, and this has the advantage that shorter 
jobs get served in less time than, say, under the first-in-first-out (FIFO) 
discipline. 

The PS discipline was introduced by Kleinrock PQ, [2], and has been the 
subject of much further investigation over the past forty years. In these mod- 
els one of the main measures of performance is a given (also called tagged) 
customer's sojourn time distribution, conditioned on the number of other 
customers in the system upon his arrival. The sojourn time is the total time 
from when a customer arrives to when that customer leaves the system, after 
being served. 

The M/M/ 1-PS queue assumes Poisson arrivals with rate A and expo- 
nential i.i.d. service times with rate fi. The traffic intensity is p = A//i. We 
shall denote the sojourn time of the tagged customer by V and the number of 
other customers present at his arrival instant by N. Then the unconditional 
sojourn time density is p(t)dt = Pr [V G (t,t + dt)'j, while the conditional 
density, conditioned on N, is p„ {t)dt = Pr [V e {t,t + dt)\lSS = n]. For the 
M/M/l-PS model we can remove the conditioning to get 

oo 

p(t) = ^(l-p)p>„(t), (1.1) 

n=0 

since N follows a geometric distribution. 

In [3] , Coffman, Muntz, and Trotter derived an expression for the Laplace 
transform of the sojourn time distribution, conditioned on both the number 
seen by an arrival and the amount of service required by the arriving cus- 
tomer, in the M/M/l-PS model. Sengupta and Jagerman [?j obtained the 
moments of the sojourn time distribution conditioned on N, and gave an 
asymptotic expansion when the number of customers in the system is large. 
Guillemin and Boyer [5] formulated Pn{t) as a spectral problem for a self- 
adjoint operator, and obtained an integral representation for the conditional 
distribution. 

Using the results in [3], Morrison f6] studied the unconditional sojourn 
time distribution p{t) in the M/M/l-PS model, in the heavy traffic limit. 
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where the Poisson arrival rate A is nearly equal to the service rate fi (thus 
p = X/fi t 1). Setting e = 1 — p, in [S] asymptotic results were obtained 
for the time scales t = 0(1), t = 0{e~^) and t = 0{e~^). Most the mass 
is concentrated in the range t = 0{e~^), and the asymptotic series involves 
modified Bessel functions. 

A service discipline seemingly unrelated to PS is random order service 
(ROS), where customers are chosen for service at random. The M/M/1- 
ROS model has been studied by many authors, see Vaulot [7j, Pollaczek [S], 
Riordan [9], Kingman [10] and Flatto [IT]. In [8] an explicit integral rep- 
resentation is derived for the generating function of the conditional waiting 
time distribution, from which the following tail behavior of the unconditional 
waiting time Wrqs is computed as 

Pr [Wros > t] ~ e-"*-^*^'^t~5/^ t oo. (1.2) 

Here a, /? and 7 are explicitly computed constants, with a = (1 — ^p)^- 
Flatto pT] obtained an integral representation for the unconditional waiting 
time distribution and derived the same tail behavior as t — >■ 00. Cohen [12] 
established the following relationship between the sojourn time in the PS 
model and the waiting time in the ROS model, 

Pr[Vps>t] = ipr[WRos>t], (1-3) 

which extends also to the more general G/M/1 case. In [T3| relations of 
the form f 1 1.3 1) are explored for other models, such as finite capacity queues, 
repairman problems, and networks. 

In this paper we study the conditional sojourn time distribution p„(t) for 
the M/M/l-PS model in two cases. First we consider a fixed p < 1 and obtain 
expansions of Pn{t) for t and/or n — 00. From these (11. 2p is readily obtained 
by using (II. ip for t large. Then we consider the heavy traffic limit where 
p t 1, and again obtain approximations for several ranges of the space-time 
plane. From these results all of the expansions in [H] can be recovered. The 
integral representation in [8] is used to derive some of the approximations. 
However, it is difficult to obtain all of the results in this paper from it. Thus, 
we derive another representation for Pn(t) using a discrete Green's function, 
which we show to be equivalent to the representation in [8]. 

We mention some related work on various PS models. In [13] we studied 
the sojourn time density conditioned on the service time in the M/M/l-PS 
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model for various asymptotic ranges, for both p < 1 and p ^ 1. The M/G/1- 
PS model was studied by Yashkov [IS], [IB], [H] and by Ott [IE]. In [IS] 
Zwart and Boxma analyze the M/G/l-PS queue with heavy tails, where the 
service density has algebraic or sub-exponential behavior. Ramaswami [20] 
studied the G/M/l-PS queue and obtained explicit results for the uncondi- 
tional moments of the sojourn time. Various asymptotic properties of the 
conditional and unconditional moments and distribution for this model were 
derived in [21]. The G/G/l-PS model has not been analyzed exactly, but 
some approximations are discussed in Sengupta [22] and the tail exponent of 
the unconditional sojourn time density was derived by Mandjes and Zwart 
[23]. A good recent survey of sojourn time asymptotics in PS queues is in 
Borst, Nunez-Queija and Zwart [23]. 

The remainder of the paper is organized as follows. In Section 2 we 
summarize and briefly discuss our main results (see Theorems 2.1-2.3). In 
Section 3 we derive the explicit formula for Pn(t) by using a discrete Green's 
function. In Section 4 we derive the asymptotic results for Pnit) for moderate 
traffic intensities p < 1. In Section 5 we consider Pn{t) for p I 1, and various 
scalings of space and time. We discuss a singular perturbation approach to 
the problem in Section 6. 
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2 Summary of results 

We consider the M/M/l-PS model with arrival rate A and we set the service 
rate = 1. Then the traffic intensity is p = A > 0. 

It was shown in [3], under the stability condition p < 1, that the recur- 
rence equation of the sojourn time density of a tagged customer, conditioned 
on the number of other customers in the system, is given by 

TL 

p'St)=PPn+l{t)-{l+p)pn{t) + ^-Pn-l{t), t>0 (2.1) 

n + 1 

with initial condition p„(0) = Taking the Laplace transform of (12.11) 

and multiplying by n + 1, we have 

{n + l)p p„+i(0) - (n + 1) (1 + p + ^) pr,{e) + n p„_i(0) = -1, (2.2) 

where pn(6') = /g°° p„(t)e"^*rft. 

Solving the recurrence equation (12. 2p . we obtain the following result. 

Theorem 2.1 The Laplace- Stieltjes transform of the conditional sojourn 
time density has the following form: 

n oo 

Pn{e)=MGnY,P'Hi + MHr, p'Gi, (2.3) 

;=o l=n+l 



where 



M = M{e)=z^(^) , (2.4) 



Gn = Gr,{9) = I z''{z+ - zY'^iz. - zr-'dz, (2.5) 







= Hn{e) = —- / z^{z+ - zy'^iz - z^r-'dz, (2.6) 

2711 

is a closed contour in the complex z-plane that encircles the segment 
Zj^] of the real axis and 

z± = z±{9) = ^\l + p + 9±^{l + p + ef-Ap\ (2.7) 

Zp L J 

a = aiO) = ^+ . (2.8) 

z<- — z_ 
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We will show in Section 3 that by taking the inverse Laplace transform, 
the conditional sojourn time density obtained from (12.31) is equivalent to the 
result in PoUaczek [8]: 



Pn{t) 



1 



27ii .A^. 



L "'0 



-(l+p-2^cos v)t jy/P^ ^) ° 



where 



(l-^e^^)™o 2^sint; 

1 



dv 



mo 



- cot V 
2 2 



(2.9) 



(2.10) 



and the contour ^* is a circle in the complex 2;-plane, centered at the origin 
and with radius less than y/p. 

Using (12. 3p and (12. 9p . we obtain the following asymptotic expansions for 
Pnit), valid for p < 1 and n and/or t — > oo. 

Theorem 2.2 For p < 1, the conditional sojourn time density has the fol- 
lowing asymptotic expansions: 

1. n ^ oo, t — s> oo with n/t > 1 — p, 

1 3p-2 



Pnit) = -A^" + - 2 

n 2(1 — pyw' 



A/-" [p(2p2 + p - 1) + 4p2Ai log(Ai) 
+ 6p(l - p)Ai - (p2 - p + 2) A?] + 0(n-3), (2.11) 



where Ai = 1 — (1 — p)t/n > 0. 
2. n ^ oo, t ^ oo with n/t = 1 - p + 0(t"^/^) 

1 



Pn(i^) 



1 — P 2-p 

n 2{i-p) 

^ V 1 +P 



X / V i-p exp , 

/o ^ 2(l + p) 



1-p 



(y-A2)'|rfi/, (2.12) 



where As = [l - (1 - p)t/n] = v^Ai = 0(1). 
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3. n ^ oo, t ^ oo with < n/t < 1 — p, 



Pn{t) ~ n-i-^ywV^ K{e,^ exp |t ( - 1 - p + + 4p) 



+ nlog 



n n 



where 



n 



0. = 9J-]=d-+Ap-l-p, 



(2.13) 



(2.14) 



a"r(a) — ^ ^-^ \ / ' — '-^ — (2.15) 
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[i-pz^Y^i + p + e 

and a = a{9) is as in l{2.8\) . 
4- n ^ oo, t — s> oo with n t^^/^ = a = 0(1), 

p,(t) ~p-"/2^(n,t) expj -(1- Vp)2t + $(n,t)} 



(2.16) 



where A{n, t) and $(n, t) have the following expressions in three ranges 
of a (as shown in Figure 1): 



(a) a > (3Vp)2/3, 
A{n,t) = 



A exp 



2(l-v^) 



( 1 - Vp) a3/4 (3 ^p{l + A) - a3/2) t 
A 



$(n,t) = (^3y^--2y/a{l + A)y'/^ 
where A = A{a) satisfies: 

I^A^I'^ a~'^l'^ = ^A{1 + A) - arcsmh(v^) 
(b) (A^/nf/' < a < (3Vp)2/3, 



(2.17) 
(2.18) 

(2.19) 



A(n,t) 



B exp 



2(1-05) 



( 1 - Vp) a3/4 ( _ 3 v/p(l - B) + a3/2) t 



(2.20) 



Mn, t) = (2-23) 



$(ri, t) = (^-3,/^^- 2v/a(l-5)) t'/^ (2.21 
where B = B{a) satisfies: 

2^5^/2^-3/2 ^ -^5(1 -5) + arcsin(\/s). (2.22 
< a < (4Vp/7r)2/3, 

(l-yp)a3/4(3v/p(l-C)+a3/2)^/^t 

$(n, t) = ( - 3Vp ^ + 2v/a(l - C)) t^/^^ (2.24) 
where C = C{a) satisfies: 

TT — arcsin(v/C') . (2.25) 

5. n = 0(1), t ^ oo, 

02/30-1/2 5/6 „-5/12 , . 

P.(t) ~ " exp { - (1 - ^p)H - 2-/^3 



;i - Vp) t5/6 

><-KT^)^i(rd)^-p(T^)^- (2-26) 



We note that for a fixed a and n, t — > oo, $(?2, t) = 0{t^^^) = 0{^/n) and 
A{n,t) = 0{t^^) = 0{n^^/'^). Despite the fact that case 4 has three different 
expressions, the functions $ and A are smooth along the transition curves 
a = nt'"^/^ = (3^)2/3 and a = [Ay/p/ir)'^/^. We also note that the contour 
integral in (12.261) is equivalent to the following infinite sum: 

-LI 1 f 1 \, _^ in + l)\ 

2m J.,, (1 - z) z-+^ ^""P U - z) ^ (/!)2n! ' 

The asymptotic sojourn time density has simpler expressions in some of 
the matching regions between the scales in Theorem 2.2. We have, between 
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cases 2 and 3, 



Pn[t) ~ — 



2ti Vl-p/V tJ \ n 

which is vahd forn — oo, t — oo, n/t < 1— p with 0(t^^/^) ^ \n/t — l+p\ <^ 
1. Between cases 3 and 4(a), we have 



2p(l-Vp)t "12(1-Vp) 

2 

X 



exp { ^ [ log(rt - 1 + 2 log(() - 3 log(n)] } , (2.28) 



which is vahd for t — > oo and 0(t^/^) ^ n ^ C'(t). Between cases 4(c) and 
5, we have 



1-7p V2, 
xexp{-(l-Vp)^t-3(|)'V^t^/^} 

^^^^{^^+2(1^}' (2.29) 

which is vahd for t ^ oc and 1 <^ n ^ 0(t^Z'^). 

By removing the condition on n, using fl2.26p in (11.11) . and noticing the 
relationship f 1 1.3 1) between processor sharing and service in random order, we 
can recover the resuhs in PoUazcek [8] and Flatto |TT], for p{t) as t oo. 

We next consider the heavy traffic case, where p is close to 1. Letting 
e = 1 — p ^ 0"^, we obtain the following results. 

Theorem 2.3 For p = 1 — e and e 0^, the conditional sojourn time 
density has the following asymptotic expansions: 
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1. n = 0(1), t = 0(1), 



Pn{t) 



-2{l-cosi;)t_ 



X ■ 



(1-e 



2 sinf 



g — TO^mo + l _j_ gTTCOtD 

where mo is defined by l\2.10\) . 
2. n = 0(1), t = a/e^ = 0{e-^), 

nVl + 



(if 



(2.30) 



p„,(t) ~ 4v^e 



3/2 _ 



exp 



3 + 4^2 



v/8 + 3a(l + 4u2) 

X exp <^ ( 7 + 3m ) cr - — 

1 /■ 1 /I 

jexp 



^8 2^"^ 2(1 + 4m2)J 



X 



27rii j^. (1 - z) 
where u = u{a) satisfies 

2u a = n lo£ 



1-z 



1 - 2iu 



dz. 



+ 



2u 



2"""Vl + 2m/ 1 + 4m2 
3. n = i/e = O(e-i), t = r/e = O(e^i), 



(2.31) 



(2.32) 



Pnit) 



T-1 At-t^ 3r2 

+ ■ 



+ 



2e 2e 



'"/« + 0(e^). (2.33) 



4. n = r]/e2 = 0{e-^), t = a/e" = 0{e~^), 

2 



Pn{t) ~ A{t], a) exp \ ^ [$(r/, f^) + | - ^ 



(2.34) 



where A{ri, a) and $(?7, a) /iaw e the following expressions in three ranges 
of the {f],cr) plane (as shown in Figure 2): 

(a) a < \ri'^l'^ - \ with r/ > 4, 



A(r7,or) 



2\lAri{\-AA) 



(1 - 4A)2 'iori\jri{\ + Ai]) - rj^ + 8ri\/r]{l + Ar]) 



X exp 



1 1. 
4 1-4A ^2 8^ . 



(2.35) 
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^{rj, a) = 3Aa - 2Jr](l + Arj) 



1-iA 



(2.36) 



where A = A{rj, a) satisfies 
2v^ 



arcsinh ( \/ Art ] + arcsinh ( 
1-AA V V 7 VV 1 - 4^ 



AA 



+ JAr){l + Ari). 



(2.37) 



(b) |7^3/2 - I < a < |7^3/2 
A{r],a) = 



f + l^-arcsin(^) 



, a>0, 



2JBr]{l + AB) 



(1 + 45)2 ^3 _ 3cTr/^r/(l - 5r/) - 8r]\/r]{l - Br)) 

7] 1 



.5 1.- 

- - ( h -)(r 

4 1 + 45 ^28^. 



X exp 



$(r/, a) = -35(7 - 2\^r]{l - Brj) + 



1 + 45 



(2.38) 
(2.39) 



where B — B{rj, a) satisfies 



2 V 5 

25^/^0- = ^ + arcsin 

1 + 45 



arcsm 



45 \ 



577(1 - Brj). 



1 + 45/ 

(2.40) 



(c) C7 > 



f + — arcsin ■ , , . , 

2 4+?7 \ y 4+r) 



2^Cri{l + AC) 



(1 + 4C)2 3(777^77(1 - C77) + 7^3 + 87;^7;(1 - C77) 

rj 1 



X exp 



- - ( h -)a 

4 1 + 4C ^28^. 



(2.41) 
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where C = 0(7], a) satisfies 
2Vc 



1 + 4C 



— arcsin I V ) arcsin 



4C 



1 + 

(2.43) 



In the heavy traffic case we can also get much more exphcit expressions 
in the matching regions. Between the regions 1 and 2 in Theorem 2.3, 



Pnit) 



2A 

Vst 
1 



exp I - 2-^/33^2/3^1/3 _1| 



X ■ 



27ri 



z) z 



n+1 



exp 



dz. 



(2.44) 



which is vahd for n = 0(1) and 1 ^ t ^ 0(e ^). In the matching region 
between cases 2 and 4(c), where 1 <^ n <^ 0{e^'^) and t = 0{e^^), 



Pn{t) 



2e2 r]-^l^u^l + Av? 
v/8 + 3^(1 + 4^2) 



exp 



(---) 



3 + 4^2 1 
2(l + 4«2) + 2 



X exp 



(2.45) 



+ 4m2 

and u = u{a) satisfies (12.321) . 

By removing the condition on n, using the results (12.331) and (12.341) in 
(11.11) . we can recover the results for p{t) in Morrison [6] for the time ranges 
t = O(e-i) and t = 0{e-^). 

Using our results for p„(t) we can also obtain some conditional limit 
laws for p{n\t) = Pn(t) (1 — p) p"'/p(t), which is the conditional probability 
of finding n other customers in the system, given the tagged customer's 
sojourn time. For t = r/e = 0(e^^) most of the mass occurs in the range 
n = ^/e = 0{e^^) and we have 

e 



p{n\t) 



(2.46) 
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where Kq{-) is the modified Bessel function. For r — oo this simphfies to 
the Gaussian 

P{n\t) ^ -=i^ exp { - - V^f}, (2.47) 



which apphes for ^ = y/r + 0(r^/''). 

When n = rj/e^ and t = a/e^, we obtain from case 4(c) in Theorem 2.3 



where C = C(cr) satisfies 



2^3/ V = TT + - 2 arcsin f W (2.49) 
1 + 4C VVl + 4C/ 

and 

g(.) ^ (l+W(l + 4a)- + 16]. (2.50) 
16[3a(l + 4C)2 + 16(l + 2C)] 

The Gaussian hmit law in (121181) apphes for r] = 4/[4C'(a) + 1] + 0{^/e). 
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3 Brief derivation of Theorem 2.1 



We use a discrete Green's function to derive (12.31) . Consider the recurrence 
equation (12. 2p . The discrete Green's function 5^(6*; n,/) satisfies 

(n + l)p$f(^; n + 1, - (n + 1)(1 + p + ^) ^(^; n, I) 

+n<^{e;n-l,l) = -5{n,l), (n, Z > 0) (3.1) 

where 5(n, /) = l[n=i} is the Kronecker delta. To construct the Green's 
function requires that we have two hnearly independent solutions to 

{n + l)pG{e-n+l,l)-{n + l){l+p + 9)G{e-n,l)+nG{e-n-l,l) = 0, (3.2) 

which is the homogeneous version of (13.11) . 
We seek solutions of (13.21) in the form 

Gn= z''g{z)dz, 

where the function g{z) and path ^ of integration in the complex z-plane 
are to be determined. Using in the above in (13. 2p and integrating by part 
yields 

Z-g(z)[pz'-{l+p + e)z+l] 

& 

- [ z"" [{pz^ -{l + p + 9)z + l)g'iz) + pzg{z)] dz = 0. (3.3) 

The first term represents contributions from the endpoints of the contour ^. 

If ( 13. 3p is to hold for all n the integrand must vanish, so that g{z) must 
satisfy the differential equation 

[pz' -{l+p + e)z + l]g\z)+ pzg{z) = 0. (3.4) 

We denote the roots of pz'^ — {l + p + 9)z + l = Oby2;+ and z^, (with 
2+ > 2_ > for real 9). These are given by (12.71) and if a is defined by (12. Sp . 
the solution for g{z) is 

g{z) = {z+-z)-''{z^-zr-\ 
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If the path of integration 3l is chosen as the segment [0, of the real 
axis, then fl3.3l) is satisfied for n > 1. Thus, we have G„ as in (12.51) . We note 
that Gn decays as n — > cxd, and is asymptotically given by 

n^oo. (3.5) 

However, G„ becomes infinite as n ^ — 1, which means that nGn-i goes to 
a nonzero limit as — 0. Thus Gn is not an acceptable solution to (13.21) at 
n = 0. 

To construct a second solution to (13. 2p . we consider another path of in- 
tegration, which is a closed contour in the complex 2;-plane, around the 
segment {z-, z+) of the real axis. Then (13. 3p is again satisfied as the endpoint 
contributions from both z = and z = z+ vanish. Thus, we have another 
solution of (13.21) . Hn, which is given by (12. 6p . if„ is finite as n —1, but 
grows as n — s> oo: 

Hnr^—-zl''-^iz+-z,r-\ n^oo. (3.6) 

Thus, the discrete Green's function can be represented by 

, . { HiGn% iin>l 
^(^;n,0 = | ^^^^^^ ifO<n<Z, (^.7) 

which has acceptable behavior both at n = and as n ^ oo. Here % 
depends only upon 6 and /. 

To determine we let n = / in (13. 2p and use the identities 

(/ + l)pGi+i = il + l)il + p + 9)Gi-l Gi^i. 

From the above we can infer a simple difference equation for the discrete 
Wronskian Gi -f^z+i — G^+i Hi, whose solution we write as 

GiHi+i-Gi^,Hi = ^^^^^, (3.8) 

where '^i = ^i(^) depends upon 9 only. Then using (13.71) in (13. 2p with n = I 
shows that % and are related by = p'~^ ^i- 
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Letting / ^ oo in fl3.8p and using the asymptotic results fl3.5p and fl3.6p . 
we determine ^\ and then obtain 

Then, we multiply (13.11) by the solution ^i{ff) to (12.20 and sum over all 
/ > 0. After some manipulation this yield 

oo 
1=0 

which is equivalent to (12.31) . 

The inverse Laplace transform gives the conditional sojourn time density 
Pn(t) as 

Pnit) = 7^ [ Pni0)e''d9, (3.9) 
where Br is a vertical contour in the complex ^-plane, with 3?(0) > — (1 — 

Now we show the equivalence between (13. 9p and (12.91) . We rewrite (12.31) 

as 

oo n 

Pn{d) = MHnJ2p'Gi + MY,p' [Gn Hi - Gl Hn] . (3. 10) 

1=0 1=0 

We deform the contour of integration in ( 13. 9p and evaluate the integrand 
along the line segments just above and just below the branch cut 3^(6') e 
[-(1 + y/py, - (1 - ^/pY]■ We denote these values of p„(6') by ^'i(6') and 
\l/2(6'), respectively. Then (13. 9p becomes 



1 /--(l-v^)' 

Pn{t) = — [vl/2(0)-x[/i(0)]e^*rf^ (3.11) 



1 f-{-^-V-p? 

and we note that ^i(^) changes to ^2(^) after making the transformation 
Zj^ — >■ Z- and a — ^ 1 — a. 

We evaluate if„ in (12. 6p by branch cut integration, which yields, for < 
a < 1, 

sin cvTT 

Hn = / r(e-^-)"-'(^+-0-"rfe (3.12) 

= z"! 2Fi{a,-n;l;- ), 

1 — a 
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where 2-^1 is the hypergeometric function. Then we observe that Hn and M 
are both invariant under the transformation z_ and a 1 — a. Thus 

we have 

^i(^) = MH^ [ ^—{z^-zy'^iz.-zr-'dz 
Jo I- pz 

+ mJ2p' iz+- zy'^iz.- z)''-^[z''Hi~ z'Hn]dz.{3.13) 
1=0 -^0 

Here we used (12.51) to evaluate (13.101) above the branch cut. But, ^'2(6') is the 
same as "^1(6), except for changing the upper hmits on both of the integrals 
in (I3.13P from z_ to z^. Also, the function {z^ — z)~°'{z^ — 2)°"^ is invariant 
under the map 2;+ — 2_, a — 1 — a. The difference is \I'2 — ^1, thus 



—iz^-zyu-zy-^dz 

pz 



+ M^^p^ Hi f ^ {z+- z)-"{z_ - zy-^dz 

- H„, I (2+ - z)-^{z_ - zr~Hz 

1 — pz 



(3.14) 



Using (13.121) and after some calculation, we find that the second part in 
(I3.14P is zero, and the integral in the first part can be evaluated by using 
contour integration (using the fact that there is a simple pole at z = 1/p). 
Thus, (13. lip becomes 

Pnit) = e''Mz-{z^-z)-'^{z-z.r-' 
2vrz Jc^ J-(i+05)2 

xT^(t^)(^) ^0^'- (3-15) 

Using the transformations Q ^ — \ — p ^ 2^/pcosf in (I3.15P and z ^ 1/z in 
(12. 6p . and changing the order of integration, we see the equivalence between 
and (12. 9p . Note that with this transformation, (a,2+,z_) becomes 
-mo,e-/yp,e--/Vp). 



17 



4 Asymptotic results for the case p < 1 



We assume that the traffic intensity p is fixed and less than one. We sketch 
the main points in deriving Theorem 2.2. We first consider n, t — >■ oo with 
n/t > 1 — p and use the result in (12.31) . From (13.51) and (13.61) . we notice that 
the first term in (12. 3p dominates the second, and thus the Laplace transform 
is asymptotically given by 

n 

Pn{9) ~ MGnY^p'Hi 
1=0 

1 ]a~l 

— — -y^^z-' (4.1) 

p{z+-z.)^ 

^ ^\''-h1^'-''Uy. (4.2) 



p[z+ -Z^) Jo 

Here we used the Euler-Maclaurin sum formula to approximate the sum in 
(14.11) by an integral. By scaling 6 = v/n = 0{l/n) and using 

^ +0{9^) as ^^0, 



1-p 



(14.21) becomes 



Pn{6) ~ -; / y^-p e ^'p^dy 

1- p ' 



= [l-{l-p)yy-''e~''ydy. (4.3) 

Then we multiply Pn{d) by e^^dO = n^^ V" an^j invert the transform to 
obtain 



1 



t 



Pn{t)-- l-(l-p)- 7>1-P. (4.4) 

nl n\ t 

To obtain the second term in (12.111) . we need the correction terms in (14. 2p . 

for which we also need the second terms in the approximations in (13. 5p and 

(IMD. 

This analysis suggests that Pn{t) is approximately zero in the range < 
n/t < 1 — p. We shall show that in this sector the density is exponentially 
small. But, we first investigate the case (1 — p)t ^ n. 
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Thus, we consider n,t ^ oo with n/t = 1 — p + 0{t ^/^). We can still use 
fl4.1l) but now scale / = y y/n = O(v^), and approximate the sum in (14.11) by 

poo 

Taking the inverse Laplace transform and scaling 9 = zu/y/n = 0{l/y/n), 
we have 

Pn{t) - — / / y^-'e-^'°^^^'^'dy)e''de 

^ n-W^ I yj^J— I e^^^'y'^dw^dy, (4.5) 



1-p 
where 



^(^'^) = 2(w)^"^ W + T^)-- 



Then by using the identity 



1 f ^Co^-+c,^dw=^^exp^ 



2vri J Br 2y/^ V 4Co 

and noting that 



n A 



2 



'n 1 — p 1 ~ P 
we exphcitly evaluate the integral over zu in (14.51) to obtain (I2.12p . 
For the case n,t oo with < n/t < 1 — p, we rewrite (12. 3p as 

OO OO 

1=0 l=n+l 

OO 



MGr^Y^p'Hi. (4.6) 



1=0 

The second sum is negligible in view of (13.51) and (13.61) . and the fact that we 
will have 6' < on this scale. The sum in (14.61) can be calculated exactly by 
using (12.61) . contour integration and the residue theorem, which yields 
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Here we used the identity pz^z^ = 1. Using fl3.5l) and fl4.7p in fl4.6p . then 
taking the inverse Laplace transform, we have 

where -F(^) = 0t + nlog(2;_(^^)). There is a saddle point at = 0* < which 
satisfies F'(6') = t + n z^(6')/2_(6') = 0, and this leads to 9^ in fl2J4D . Hence 
using the saddle point method leads to (12 ■131) . 

The expression f l2.27p in the matching region between cases 2 and 3 follows 
by letting A2 —00 in (12.121) . or letting n/t ^ 1 — p in (12.131) (which 
corresponds to 6'=,, ^ 0). 

Now we consider n, t — > 00 with nt^"^^^ = a = 0(1). We first note that 
0^ ^ 9^ = -{1 - as njt in (12^4]) . Expressions (gl]) and (gTD are 
still valid and we have 

v — p z_\z_ — \) 

Then, by taking the inverse Laplace transform, the conditional sojourn time 
density is asymptotically given by the double integral 

p„W~^/ ^(^X { r (4.8) 



2Tii J^^. 1 ~ p z„ \z_ — 1 
Scaling 9 = 9c + s/n (with s > 0), we notice that 



z± = p '/'±p-'/'^ + 0(-). 



Thus, we scale z = z^— p '^^^ y/ ^/n (with ?/ > 0). The inner integral in (14. 8p . 
which is G„, is asymptotically equal to 



Gn- r^=pL=exp{V^c^{y,s) + -^- ^y^'^'^' ]dy, (4.9) 
Jo vW + 2V^ 2^ 2^ J 



where 

,1/4 



2vs \y + 2^ysJ 
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The function (j){y, s) has its maximum at y = y*{s) = ^/s'^T'^ — a/s > 0, 
which satisfies 4>y{y*,s) = and (pyyiy^.s) < 0. Hence, using the Laplace 
method in fl4.9p . we have 



[n {s + 7p)] 
In the hmit of 6 ^ 6c we have 



1- pz^\z^ - iJ yp(i-^) VI -^y 

with which Pnit) in (14. 8p asymptotically becomes 



Pn(0 ~ T^-^- — exp <^ — I p-i^-Vp) * 



X 

27ri 



— / (s + v/;^)-^/^ev^*°(^)ds, (4.11) 



where $o(s) = 0(?/*(s),s) + a~'^/^s, and Br' is a vertical contour in the 
complex s-plane. The saddle point s=k = s*(a) satisfies $o(s*) = 0, which 
implies that 



4pV4,-3/. ,3/2_^ (^^^^-^) -27ZMV^ = 0. (4.12) 

If let s=K(a) = y/pA > 0, fl4.12p is equivalent to (12.191) . Using the saddle point 
method in (14.111) . we obtain (12.161) with A{n,t) and $(n, t) = ^/n(^o{s^) as 
in (EHZD and fl27[8D . We note that a ^ (Syp)^/^ as A and that A is an 
increasing function of a, so the above result is valid for a > (Sy/pY^^. 

Alternately, on the scale n = 0{t^/^) we use the representation (12.90 with 
the scaling z = y/p{l — t~^^^w) and v = t~^^^u {w > 0,u > 0). Then we 
have 

z-"~p-"/2exp (^awt'/' + ^), 

[1 + p- 27p cos(i;)] t ~ (1 - 7p)2 t + y^u^ t^/\ 
sin(w) / TTt^/^- 



t"^/^ exp — y 



]^ _^ g7rcot{i,) U 
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: exp ' 



and 



It follows that 

2p-"/2 / 



x( [ g{u,w)e-''^'f^''^'"Uu)dw, (4.13) 



where 



1/ r ~i 

q(u, w) = = exp < — — >, 

^ ^ + m2 ^12(m;2 + m2)/' 

2 , ^ ^ T fw -iu 

f[u,w) = ^u + — - 

■u 2^ V-i 



■u 2u \w + 

and the contour Br^, is a vertical contour in the complex w-plane with 
sufficiently large. The function f{u,w) has its maximum at = u^:{w), 
which satisfies 

d s ^ TT i ^ fw — iu\ w 

—f{u, w) = 2^pu - - + — log —— - = 0. 4.14 

ou Iw^ \w + iuJ u^w^ + w^) 

Then using the Laplace method in the inner integral of fl4.13p implies that 



Vn{t) ~ -^p-"/2exp(-^)r^/ee-(^-v^)^* 

1 - ^^0 VI - ./pj 



271"^ J Br, \/ fuu{.U^,w) 



"^kw. (4.15) 



Let h{w) = aw — /(n*(w),ty). Then the saddle point equation h'{w) = 
along with fl4.14p leads to 



wl + ul = -. (4.16) 
a 
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Here we denote the solution to h'{w) = by = Wo{a) and set Uq = u^:{wo). 
Applying the saddle point method to fl4.15p . we find that 



(1 - VP) V fuu{uQ,WQ) h"{wo) 
If we let C = 1 - awg (0 < C < 1), then from (14161) it follows that 

Wo = V{l~C)/a > and Uq = ^/Cj^ > 0. (4.18) 

Since Wq and Mq satisfy fl4.14p . using (14.181) in (I4.14p leads to (I2.25p . Hence, 
after some simplification in (14.171) . we obtain (12.161) with A and $ = t^^^ h{wo) 
as in fl^:^ and (^M>- We note that a^OasC^O, (A^/n)^/^ as 
C 1 and C = C(a) is an increasing function of a. Thus the above result 
is valid for < a < (4^/7?)^/^. 

We now consider the range (4y/p/7r)^/^ < a < (3y^)^/^. This is difficult 
to treat using either of the representations in ( 12. Sp and ( 12. 9p . as the various 
saddle points become complex. However, we now show that the results in 
case 4(b) of Theorem 2.2 can be obtained by smoothly continuing the results 
for case 4(a), or those of case 4(c). 

First we consider (I2.19p which we rewrite as 

„ ,„ , arcsinh (v^) 

2VpAa~^^^ = ViTA (4.19) 

V A 

= ^A + 0{A'), 

where the right side of ( 14.191) is an analytic function of A. The curve a = 
(3yp)2/3 corresponds to A = 0. Setting A = -B in (ICTD we obtain (12:221) . 
which is the analytic continuation of (I4.19P into the range A < 0. Then (I2.20p 
and (12^211) follow by replacing A by -5 in (12X71) and ([238]). We now show 
that case 4(b) also follows by the continuation of case 4(c), as a increases 
past (4y/p/7r)^/^, which corresponds to C = 1. The smooth continuation of 
( I2.25P as C increases past C = 1 follows by replacing C by i?, aTcsm{\/C) 
by TT — arcsin(\/5) and a/1 — C by — Vl ~ B. Note that viewing C as a 
function of {n, t), 1 — C has a double zero along nt'"^^^ = (4^/p/7r)^/^. These 
observations show that the three cases in Theorem 2.2 for n = 0{t^^^) really 
correspond to a single asymptotic scale. A geometric interpretation of these 
three cases is given in Section 6. 
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In the matching region between cases 3 and 4(a) in Theorem 2.2, we let 
n/t — > in (12.131) . which yields (12.281) . On the other hand, letting a — oo in 
(lilHD, we obtain 

A=^a'-3 log(a) + (l + log(p)) + O (i) . (4.20) 

Using f|4:2nil in fl2Tfl) and fl2J8D also leads to (Km . 

Finally, we consider n = 0(1) and t ^ oo. We use the representation 
( 12:91) and scale v = t'^/^ V. Then the inner integral in (12. 9p becomes 



(1- Vp)(v^-^) 

Using the Laplace method we find that the integrand is maximal at = 

1 /3 

(27^) ' ^^"-^ then making the transformation z in the outer integral 

in ([2l9D lead to fl2:26|) . 

To verify the asymptotic matching between cases 4(c) and 5 in Theorem 
2.2, we let a ^ in (l2^ . It follows that 

C~(— ) a. (4.21) 

Using flT2TD in (12:231) and KJIM) yields fl2:29|) . On the other hand, we can let 
n — > cxD in (I2.26P . We scale z = 1 — w/ ^Jn in the contour integral in (12.261) 
and use the saddle point method. There is a saddle point at w = 1 and we 
obtain 



exp 



27ri 7^. (1 - z) \\-z) l^n^l^ 

This also leads to (I2.29p . which verifies the matching. 



dz ^fi-— e^v^. (4.22) 
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5 Asymptotic results for the case p ~ 1 



Now we consider the case in which the traffic intensity is close to one. Letting 
p = 1 — e with < e -C 1, we sketch the main points in deriving Theorem 
2.3. 

First, we consider n = 0(1) and t = 0(1). Replacing p by 1 in (12.91) 
leads to (12.301) . On this scale the solution does not simplify much, but there 
is little probability mass in heavy traffic on the time scale t = 0(1). 

Next, we consider n = 0(1) but very large time scales t = a/e^ = 0(e~^). 
We use (12.91) and scale v = 0(e). Then (12.311) is obtained by making the 
transformation z —>■ ^/pz in the outer integral and using the Laplace method 
in the inner integral, where the major contribution comes from the point 
u{a), which satisfies (12.321) . 

To verify the matching between cases 1 and 2 in Theorem 2.3, we let 
cT ^ in (E32D, which yields 

Using (15. ip in (I2.3ip leads to (I2.44p . On the other hand, we can let t ^ oo 
and scale v = 0{t-^^^) in flOOD . Then using the Laplace method in the inner 
integral also yields (12.440 . This implies that there are no other time scales 
between t = 0(1) and t = 0(6"^). 

For the case n = ^/e = 0(e^^) and t = r/e = 0(e~^), by similar argu- 
ments as in the case when p < 1 and n = 0{t) in Section 4, we can easily 
obtain (12331) . We note that if ^ ^ cx) and r ^ oo but ^/r = 0(1), fl233l) is 
still valid. 

For t = 0(e~^), we rewrite (I2.33P as 

p„(t) = e Po(e, r) + P^{^, r) + 0{e^), (5.2) 

and remove the condition on n by using (15.21) in (II. ip with the scaling n = ^/e. 
It follows that 



e-« Po(e, r)d^ + e'^ [p^i^, r) - | Po(e, r) 



2eKo(2v^) + - (6-r)iro(2v^)- v^Ki(2v/7) 
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Here Kq and Ki are the modified Bessel functions. Tliis recovers tlie result 
in Morrison [B] for tlie time range t = 0(e~^), after we take into account that 
the results in [6] are for Pr [V > t] . 

Now we consider larger space-time scales, with n = rj/e^ = 0{e~^) and 
t = a/e^ = 0{e~^). We use a similar method as for p < 1 with the scale n = 
0(t^/'^). However, the heavy traffic assumption changes some of the saddle 
point calculations. We note that fl4.8p is still valid in the heavy traffic case, 
but 9^ = -eV4 + 0(e3) as e ^ 0. Then we scale 6 = (3 (-1/4 < /3 < 0) 
and notice that 

z± = l + e + 0(e ). 

Thus, by scaling z = 1 + re, the inner integral in (14. 8p . which is is 
asymptotically given by 

wi+s 

Gn~ / ' Jo(r,/?)e'^(^'^)/^tir, (5.3) 



where 



, ^ 2/3+1 

1 - 2r - ^Jl+Afj\ 2{4;3+l)3/2 



v/(l -2r)2 - (1 + 4/?) V 1 - 2r + ^/TTW 
r (1 + 3/3)+/? (1 + 2/3) r^r^ 



X exp 



;i + 4/3)(r2-r-/3) 



ip[r, p) = rr] -\ log 



and 

il](r fi^ = r n -i- 

The major contribution to the integral in (15.31) comes from r^, = r^{(3), which 
satisfies ipr{i^*, P) = 0, so that — r^: — {(3 + 1/r]) = 0. For ip to be maximal 
we need 

^rr{r*,f3) = rf (2r, - 1) < 0, 

and this implies that 



1 - Vl + 4/3 + 4/r/] < ^. (5.4) 



Using the standard Laplace method in (15. 3p . we obtain 

V2^Jl(/?)e'^(^••'^)/^ (5.5) 
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where 

Jo(r.(^),/3) 



r/ (1 + 4/3 + 4/r7)V4- 
Using (15. 5p in (I4.8p . and the expansion 

, , 1+2/3 

1\° 1 f 1 + + ^VT+w 2(1+4/3)3/^ /1 + 2/5 



we have 



27ri 

where 



o 3/2 /• 

i'n(^) ~ / -^il/^) '^2(/5) e*('^)Ad/3, (5.6) 

/Br" 



1+2/3 

1 /l+V^rT4^\ 2(1+4/3)3/2 /I + 2/3 

•^2(/3) = ^= ^ , exp 



^\\l-^JlT^J ^Vl + 4/5 

^(/3) =/3a + -=i= log n^'^j-^ l +^(r,(/3),/3), 



v/rT4;5 ^Vi-v^rT4;5 

and -Br" is a vertical contour in the complex /3-plane. Then \l/'(/3) = implies 
that there is a saddle point at /3o = f^oij], a), which satisfies 



yi + 4/?o 

2/3o(l - 2ro) 



2ro + a/3o(l + 4/?o)(l - 2ro) - (1 + 4/?o)(?//5o + 1) 



- log + log f l±v™) . (5.7) 

^ Vi-2ro + v^r+4:3^; ^ Vi- v^r+4:3^/ 

Here we denote r^,(/5o) by tq. Thus, from (15.61) . the saddle point method 
imphes that 

V^"(/3o) 

If we let /?o = ^ — 1/4 (0 < A < 1/4), then after some simplification, we have 

HMM^^l^,,,) (5.8) 



and 



vi/(/5o)=<f(r7,a) + |-^, (5.9) 
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where A and $ are given by fl2.35p and 02.361) . Equation fl2.37p is derived by 
using ([53D and ([52D- We note from ^M) that a T - | as 1 i 0. This 
imphes that rj > 4, since a > 0. Hence (12.351) and fl2.36p are only vahd in 
the range a < |?7^/^ — | with 77 > 4. 

Ahernately, on the scale n = 0(e^^) and t = 0{e^^), we use the repre- 
sentation in (12. 9p and scale v = and 2; = ! — ea;(7>0,a;>0). Then we 
have 

(1 - ^e'")"'" 2 f ^ 1 / l-2z7 \ 47^ + 3 ) 

(1 - ^e-™)'"«+i ^ ev/472 + 1 I 2e7 °^Vl + 2i7v' 2(472 + 1)/' 



-2 



mo+l 



exp 



^472 + {2uj - 1)2 

472 - 4w + 3 



l2e7 



2cj - 1 



2^7 \ 



2cu - 1 + 2^7/ 



X exp 



2[472 + (2cu- 1)2] 
and (12. 9p is asymptotically given by 



where 
/(7,w) 



^y^^exp(!^ + !^)( / /(7,..)e-^(--)Arf7)rfc., (5.10) 



(472 + 1) [472 + {2u - 1)2] 

472 - 4cj + 3 472 + 3 



A;(7,a;) 



X exp 
1 



2[472 + (2cu- 1)2] 2(472 + 1; 

2cj - 1 - 2i7\ 



r 
2 



O" 



TT 



7 
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log 



2cu - 1 + 2^7/ 



+ log 



1 - 2z7 
l + 2i7/J 



and Br is a vertical contour in the complex tu-plane with ^{u) sufficiently 
large. Thus, the major contribution to the inner integral in (15.100 comes 
from 7* = 7*(co'), which satisfies ^^(7^,0;) = 0, that is 



2^7^ - 



4cu7, (472 + 2CJ-I) 
(472 + 1) [472 + (2a; -1)2] 



TT 



+2 



log 



2a; - 1 - 2z7, 
2c; - 1 + 2^7, 
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log 



1 - 2z7, 
1 + 2z7, 



(5.11) 



Using the Laplace method in the inner integral in flS.lOp yields 

pM^^^I e-'-'/'e^-'-Vrfv. (5.12) 

^TT^ J Br ^/k^^{'y^,uj) 

where = rju — k{j^:{uj),uj). The saddle point equation ^[{uj) = has 

the solution luq = ujo^r], a), which satisfies 

47o' + (2u;o - 1)' = ^, 7o = 7*M- (5.13) 
In order that ^i^co-o) > 0, from 05.131) it follows that 



^o = l + ^l^-ll (5.14) 



Then using the saddle point method in (I5.12p . we have 



(5.15) 

If let 7o = VC (0 < C < 1/r/), then from flCTj) we obtain cjq = 1/2 + 



\|^-C. Using this in flSTTTj) leads to ([233]). It follows that 

A(r7,a) = ^ =6^^-°/^ 
A/A;T,T,(7o,t^o) '^''/(t^o) 

and 

$(r/,a) = $i(^o)-| + ^ 
in flCTll and (ElSl). We note from (ICTD that at (5 = I/77, we have 
1 



vr 4^ . ( / 4 

— I arcsm 



2 4 + ^ VV4 + r7 



(5.16) 



We also have a t^/{2C^/^) as C ^ 0. Hence (ICTD and (12:^21) are only 
valid when a exceeds the right side of (15.161) . 

The range of a between 

- = \^"'-l (5-17) 
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and fl5.16l) is difficult to treat using either of the representations in fl2.3p and 
(12.91) . as the various saddle points become complex. Similarly as in case 4(b) 
in Theorem 2.2, we now show that the results in case 4(b) of Theorem 2.3 
can be obtained by smoothly continuing the results for case 4(a), or those of 
case 4(c). 

First we consider (I2.37p which we rewrite as 



oT 2 1 / 1 / I AA 

2 A a = ^ arcsinh lyAr]] H t= arcsinh I i 



1 - AA ^yx ^ ' yA V V 1 - 



+ ^r]{l + Ari) (5.18) 

'2 n/9 16\ 



{{^"'"-i)^ + 0{A') 



where the right side of flS.lSp is an analytic function of A. The curve (15.170 
corresponds to A = 0. Setting A = —B in (15.181) we obtain (12.401) . which is 
the analytic continuation of (15.181) into the range A < 0. Then (12.231) and 
flT^ follow by replacing A by -B in ([235D and flT^ . We now show that 
case 4(b) also follows by the continuation of case 4(c), as a decreases past 
the curve (I5.16p . which corresponds to C = l/rj. The smooth continuation of 

(I2.43P as C increases past C = 1/r] follows by replacing C by B, arcsin(-y/ C rf) 



by 7r — arcsin( Y B rj) and yl — Cr] by —y 1 — B rj. Note that viewing C 

as a function of {ri,a), 1 — C t] has a double zero along the curve (I5.16p . 
These observations show that the three cases in item 4 of Theorem 2.3 really 
correspond to a single asymptotic scale. A geometric interpretation of these 
three cases is also given in Section 6. 

Now we consider the matching between cases 3 and 4(a) in Theorem 2.3. 
If we fix a but let rj = (/e ^ oo in (12.371) , it follows that 

A=\-'-e-^^e' + 0{e% (5.19) 

Using flCT]) in fl2:35|) and fl236D . we have 
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Alternately, if we fix rj but let cr = re^ — > in fl2.37p . we have 
Then p„ (t) becomes 

V ^ 
These calculations verify the matching. 

Next, we consider the matching between cases 2 and 4(c) in Theorem 2.3. 
Since 70 and Uq satisfy flS.llI) . we use (15.141) in (15. lip and let 77 — 0. Then 
the leading term in the asymptotic expansion in (15. lip leads to (I2.32p with 
u = 7o. Then letting 17 in ([MID and (12:421) leads to ([235]). On the other 
hand, if we scale z = 1 — ew in (I2.3ip and use the saddle point method, we 
obtain 

Using (I5.2ip in (I2.3ip also leads to (12.451) . This verifies the matching. 

To remove the condition on n, we use (15.151) in (11.11) . Since p" ~ exp(?7/e— 
77/2), it follows that 

r ''-'"■"t, , exp{-ln(,..)+(M-|)}d,. (5.22) 

Jo ^/k^^{-fo,uJo)^'i{uJo) l- e V 2 2/ J 

Here fl{ri,a) = r] — $i(u;o(r7, a)). Then ^^(r], a) = implies that the major 
contribution comes from 77 = r]Q{a) which satisfies uJo{r]Q, a) = 1. Then using 
the Laplace method in (I5.22p . we have 



^3/2 V^TT ^(70;^0 



p{t) ~ e^/^ "^'"'^"^ e-'"^vo,^)/e_ (5 23) 



Here we set Uq = cuq (770,0") and 70 = %{ujQ{riQ,a)). From (15.140 . if we let 
7o = I cot (^) {0 < t/j < it), then r^o = 4sin^ (^). Thus, (15. lip becomes 

cr = or(^) = 4( sin(?/') + ij) tan^ (^^ j 

and we also have 

n{r]o,cT) = -2^ tan (^) + ^ csc^ = Fo(^). 
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Letting 




(I5.23P becomes 

which recovers the result in Morrison [6] on the time scale t = 0(e~^). 
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6 Singular perturbation method 

Now we discuss an alternate, singular perturbation approach for deriving the 
asymptotic approximations. 

We first assume that the traffic intensity p is fixed and less than one. We 
introduce a small parameter 5 (0 < 5 ^ 1), let n = N/5, t = T/6, and 
expand Pn(t) as follows: 

Pr.it) = 6qo{N, T) + 6'q,{N, T) + 0{6'). (6.1) 

Using the recurrence equation fl2.ll) . the leading term in fl6.ip satisfies 

with the initial condition qQ{N,0) = 1/N. Solving this PDE by the method 
of characteristics, we obtain 



*(iV.T) = l[l-(l-p)l 



1-P 



(6.2) 



Using (16. 2p . we can solve for qi = qi{N, T), which satisfies the PDE 

dqi _ 1 + pd^qo 1 dqo ]^ < J_ _ 1_ 

dT ~ 2 dm NdN ^ ^'dN N^"^"" AT^^' 

with the initial conditon gi(iV, 0) = —l/N"^. Hence, we have 



qi{N,T) 



1 



3p-2 



A,'-" p{2p' + p-l)+ 4p'A, log(Ai 



2(1 -p)3iV2 

+ 6p(l-p)Ai-(p2_p + 2)A2 



(6.3) 



Using (16.21) and (16.31) in (16. ip . we get (12. lip upon setting 5 = 1, so that 
{N,T) = {n,t). From (16. 2p we note that this result is valid for n, t ^ oo 
with n/t > 1 ~ p. 

Next, we consider n, t — > oo with < n/t < 1 — p. We assume that Pn{t) 
has an expansion in the following form 
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Using this in (^Jij yields the PDEs 
and 



(6.4) 



K 



(1) 



'(1) 



+ (f e^v^ + i e^viv e-«- - ^ e-«-) iT^^). (6.5) 



The PDE (16. 4p can be solved by the method of characteristics with all of the 
rays starting from the origin {N,T) = (0, 0). This leads to 



eiN,T) = T{-l-p+\l ^ + 4p]+N\o^ 



1 / N m , 



. (6.6) 



The function K^-^'' (N, T) cannot be determined completely, but from (16.51) 
we find that it has the form 



(6.7) 



Since Pn{i) must ultimately be independent of 5, we can set 5 = 1. Thus, we 
p„(t) ~ n-i-iVi+4p*Vn^ k{^) e®("'*), < ^ < 1 - p. (6.8) 

On the scale n,t —>■ oo with n/t = 1 — p + 0{t^^^'^), we let = (1 — 
p)T + V6S. From (ES]), it follows that 



1 

N 



Then we expand Pn{t) in the form 

p„(t) ~ A^"^^(A^, 5). 
Using (16. 9p in (12. ip . we obtain for the heat equation 

^^5S = (1-P)^7V. 



(6.9) 



(6.10) 
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We next obtain two matching conditions between the scale S = 0{1) and the 
ranges n/t > 1- p and < n/t < 1 - p. We first let = (1 - p)T + S 
in (16.11) and fl6.2p . which yields 



2-P 



If (16.11) and (16. 9p were to match, the above should agree with the behavior 
of (16. 9p as 5 — oo, which implies that = 2(1- p) ^^"^ 



.^{N, S) ~ S~ as S ^00. 



(6.11) 



We next consider the matching between the scales n/t ^ 1 — p and < 
n/t<l- p. We let n = N/6 and t = T/6 in 1^ and let N/T ^ 1 - p. In 
this limit we have 



e(,M)4ewT)^-i|j^ 



1 2 



so that 



Pn{t) 



3-p 
2{l-p) 



exp 



2{1 + p)N 



K 



N/T^l-p 



We furthermore assume that K{-) has some algebraic behavior as N/T 
1 — p, in the form 



K 



Co 

Co 



(1-p)- 



N' 



T 

l-p)V6 



S_ 

iVJ 



N 

as — 
T 



where Cq and vi are constants that will be determined later. We thus obtain 
the second matching condition 



Co 



3-p 
2{l-p) 



P 



S_ 

N 



^1 



exp 



(I-P)g^ 
2{l +p)N 



, (6.12) 



as 5" — >• —00. By comparing powers of 5 in (I6.12p . it follows that 

3-p vi 
- 2(W) 2" 
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and thus vi = The matching conditions also suggest that we seek a 

solution of fl6.10p in terms of the similarity variable S/ VN, with 



^{N, S) = N''' -^{-1=) = ^''^ ^(^2), (6.13) 



S_ 

. ~N 
where 

1 3 — p z/i p 

^ w ~ 2(1 -p) " y = 2(1 -p)- 

Then (16.111) and f l6.12p give the behavior of .^(A2) as A2 — > ±00, as 

^(Aa) ~ A2^ as A2 ^ +00, 

jr(A2) ~ Co(l -p)"T^ (-Aa)"^ 

X exp I - ^l^A^} as A2 ^ -00. (6.14) 

Using fl6.13p in (16.100 . we find that satisfies the parabolic cylinder equa- 
tion 

(1 + p) ^"(Aa) + (1 - p) A2 ^'(A2) - p ^(Aa) = 0. 
The solution that satisfies both matching conditions is given by 

Combining (16.90 . (16.130 and (16.150 and setting 5 = 1 we regain (12.120 . We 
also note that by letting A2 — > —00 in (16.150 and using (I6.14p . we determine 
Co as 



27r V 1 + P VI -p 

Now we consider n,t 00 and nt~'^^^ = a = 0(1). From the case 
n,t —>■ 00 with < n/t < 1 — p, if we let ra/t ^ in (16. 6p . it follows that 

^ein,t) ^ p~n/2 I _ (1 _ ^ft + Oin'/t)^. 

Then we set 

Pn(t)=p-"/'e-(i-v^)'*i?„(t). (6.16) 
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Using fl6.16p in fl2.ll) . we have the following recurrence equation 

1 Tl 

—Kit) = -2 Kit) + ^^i?„_i(t). (6.17) 

^/p n + 1 

Letting n = Y 5'"^^^ and t = T/5 we assume that Rnit) has the following 
asymptotic expansion: 

Kit) ^ 6"' AiY,T) exp{r^/3^(r,T)}. (6.18) 
Using (I6.18P in (I6.17p . the perturbation method yields the PDEs 

^vl/^ = (vI/^)2 _ 1 (6.19) 

and 

At - 2^y Ay = (^^yy + — ^y) A. (6.20) 

We use the method of characteristics to solve fl6.19p . with all of the rays 
coming from the origin iY,T) = (0,0). We find that the geometry of the 
rays naturally defines three regions in the (y, T) plane (as shown in Figure 
1). The first region corresponds to \I/t > and \E'y < 0, where the rays A> 
satisfy fl2.19p . Here \E't = ^/pA/Y is constant along a ray. When A = 0, 
a = nt~'^^^ = (3^/p)^/'^, which is the ray denoted by the dashed curve in 
Figure 1. The second region corresponds to \1/t < and \l/y < 0, where the 
rays < 5 < 1 satisfy f l2.22p . with now \1/t = —^B/Y being constant along 
a ray. If let 5 ^ 1 in ([222), then a iA^/n)^/^. a = i^^jufl^ is not a 
ray, which we denote as a dotted curve in Figure 1. This curve corresponds 
to the locus of the maximum values of Y achieved along the rays that start 
from (0, 0) and return to F = at some later T > 0. The third region 
corresponds to \I/t < and \E'y > 0, where the rays < C < 1 satisfy fl2.25p . 
Thus, (12.181) . fl2.2ip and (12.241) are obtained from the corresponding region 
with 5-i/3^(r,T) = $(n,t). 

The function A(y, T) cannot be determined completely from fl6.20p . but 
we find that it must have the form 

A(r,T) = lAo(^)=^Ao(a), 

We shall obtain the behaviors of Ao(a) as a ^ and a — > oo below. Thus, 
by setting 5 = 1, p„(t) has the following asymptotic approximation: 

Pn(t)~^Ao(a) expj -(1- V;^)2t + $(n,t)|, (6.21) 
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where $(n, t) has three different forms, but passes smoothly through the two 
transition curves in Figure 1. 

For the scale n = 0(1) and t — oo, we assume that Pn{t) has the form 
Using the above expansion in fl2.1l) . we obtain the difference equation 

- + ^T^n-X = 0. (6.22) 

n + 1 

Using generating functions we can express 0^n{T) in terms of I^q{T), as 



1 C e 



-1 



1 



^„(T) = ^o(T) <p —-n re—dz. (6.23) 

Then, by setting 5 = 1, we have for n = 0(1) 

1 f e-i 



27r7 Jc^, (1 - -2) 



In (16.81) . (16.211) and (16.241) . the functions K, Aq and cannot be com- 
pletely determined by the perturbation method. But by examing the match- 
ing between the scales, we can find their structure in the matching regions. 

We first consider the matching region between the scales < n/t < 1 — p 
and n = 0{t^/'^). If let n/t ^ in ([6SD, we have 

^-i-iVi+4ptVn%e(n,t) _ P:!^g-(i-v/p)^t expj _ ^ _ /pliog(n)| 

n I n ) 

and we assume that K{n/t) has the form 

If let a = nt^"^^'^ — ^> oo in (I6.2ip in the first region where a > (3y^)^/^, we 
can use (I4.20p and obtain 

$(n, t) ~ -^-^ + v^- [ log(p) -1 + 2 log(t) - 3 log(r7 
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Then we assume that Ao(ci) has an algebraic behavior as a ^ oo, in the form 
Ao(aj ~ C2(^^^j as a = nt ' oo. 

If these two scales are to match in an intermediate limit where n/t ^ and 
a — i> oo, it follows that 

Thus, we conclude that Ci = C2, /3i = 1 and (32 = 0. We also conclude that 
71 (n, t) = [ log(p) - 1 + 2 log(t) - 2 log(n)] . 

We now consider the matching region between the scales n = 0(t^/^) and 
n = 0(1), t — s> 00. If we let a — s> in (I6.2ip in the third region where 
< a < (4y/p/7r)^/^, we have fl4.2ip . which implies that 



$(n,t)~-3(^) pi/6ti/3 + 2v/^. 

We then assume that Ao(a) has the form 

Ao{a) c-if^^Y' as a = nr^/^ ^ 0. 

We let n — s> 00 in (16.241) and use (14.221) . and also assume the following form 
for ^o{t) 

^o{t) ~ C4 1^* exp (74 1^/^) as t ^ 00. 

Note that the exponential factor is indicated by the behavior of $ as a ^ 0. 
Then the matching holds provided that 

Thus, we conclude that 

/33 = -1/4, = -5/6 and 



74 = 
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K2) " 



We have thus shown that much, but certainly not all, of the asymptotic 
structure of p„(t) in Theorem 2.2 can be obtained by perturbation methods, 
which make no recourse to the exact solutions in Theorem 2.1 and (12. 9p . We 
can use this method to obtain the full asymptotic series for the scales for 
cases 1 and 2 in Theorem 2.2. For the scales in cases 3, 4 and 5 we can 
obtain partial information only. Specifically, we can get (12.131) only up to the 
unknown function K{n/t), for which we can infer the behavior asn/t 1— p 
and as n/t ^ (up to the constant ci). Of course this function was fully 
determined by the saddle point method, as given in (I2.15p . From (12.151) we 
can easily show that the behaviors as n/t ^ 1 — p and as n/t — > were 
correctly predicted by the matching arguments. 

Similarly, the perturbation method yielded the function ^{n,t) in (12.161) 
completely, for all 3 ranges of a. But, only partial information could be 
obtained about A{n,t) = t~^Ao(a). Specifically, we obtained the behavior 
of Ao as a — > oo (up to the constant C2) and as a — >• (up to the constant 
C3). For t ^ 00 with n = 0(1) we could determine the expansion of Pn{t) 
up to the constant C4, which is expressible in terms of C3. For the scale 
n = 0{t^/^) the perturbation method led to a nice geometric interpretation 
of the 3 sub-cases in item 4 of the Theorem 2.2. 

Next we very briefly discuss the heavy traffic case via perturbation ex- 
pansions. We set p = 1 — e with < e -C 1, and first consider n = ^/e and 
t = r/e, expanding Pnif) as follows: 

Vnit) = ePo{i,T) + Pr{i,T) + 0{e'). 

Using (I6.25P in (12.11) we obtain to leading order 



(6.25) 



dr 



with Po(e,0) = - 



so that 



To obtain Pi(^, r), we need to solve 



(6.26) 



(6.27) 



dPi d^Po 



1-1 



with the initial condition Pi(^,0) = —1/^^, and hence 

3r2 ^ r3_ 



Pi 



A(e,r) 



1 4r - 
- + 



2^3 



2e 



(6.28) 
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Using (Km and (KM in fl^:^ leads to ([233]) • This result is valid for 
n, t — ^ oo with n/t = ^/r = 0(1). 

Next, we consider the scale n = rj/e^ = 0(e^^) and t = aje" = 0(e^^). 
We note that from (16.211) . by setting p = 1 — e, we have 

p-"/^e~(-v^)^* = exp{l(|-9+0(l)}. 
This leads us to seek the asymptotic expression of p„(t) in the form 

Pn{t) ~ e^^ A{V, ^) exp |i [$(r/, a) + | - ^] }. (6.29) 
Using firai) in (O), we have the following PDE: 

= -. (6.30) 

This PDE is essentially the same as that in (16.191) . However, now we impose 
the initial condition $(77, 0) = —ri/2, which is necessary if the expansion in 
(I6.29P is to match to (16.271) as r/^ = a/(^er]) — >• 00. Thus we must solve 
(I6.30p using characteristic curves (rays) that start from (?], a) = (?7,0), with 
rj > 0. We again find that the geometry of the rays naturally divides the 
(?7, cr) plane into 3 parts (see Figure 2). In Figure 2 the rays in region (1) 
always have da/drj > and these never hit 77 = (the scaled time axis). 
Regions (2) and (3) are filled by rays that do hit t] = 0, and region (2) has 
da/drj > along a ray, while region (3) has da/drj < 0. 

The dashed curve is (I5.17p . which is a ray corresponding to $0- = ^ = 
in (12.371) . that separates regions (1) and (2). Letting B 1/rj in (12.401) . we 
obtain (15.161) . which is not a ray, and is shown by the dotted curve, which 
also separates regions (2) and (3). Expressions (12.361) . (12.391) and (12.421) 
are obtained in the corresponding regions by solving (I6.30p . The function 
A(?7, 0") cannot be determined completely by the perturbation method, but 
some partial results can be obtained by using matching arguments, as was 
the case when p < 1. 
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Figure 1: The rays in the (Y, T) plane for the p < 1 case. 
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Figure 2: The rays in the {rj, a) plane for the heavy traffic case. 
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